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We investigate the effects of pliase-breaking events on electronic transport through ballistic 
chaotic cavities. We simulate phase-breaking by a fictitious lead connecting the cavity to a phase- 
randomizing reservoir and introduce a statistical description for the total scattering matrix, includ- 
ing the additional lead. For strong phase-breaking, the average and variance of the conductance 
are calculated analytically. Combining these results with those in the absence of phase-breaking, 
we propose an interpolation formula, show that it is an excellent description of random-matrix 
numerical calculations, and obtain good agreement with several recent experiments. 
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Recently there has been great interest in the effects of 
quantum-mechanical interferencje on electronic transport 
through ballistic quantum dots.EJ In these microstructures 
both the phase-coherence length and the elastic mean free 
path exceed the system dimensions. Thus the leads into 
the dot can be thought of as electron waveguides and the 
dot itself as a resajiaiji|t.j2aYitir. 

ExperimentallyoBu'ErElu'ErEl one observes random (but 
reproduciblft^ iluctuations in the conductance as the mag- 
netic |field,BQQ the Fermi energy,El or the shape of the 
cavityQ is changed. The sensitivity to small changes in 
these parameters shows that these fluctuations are caused 
by quantum interference. In addition to conductance 
fluctuations, several interference effects which surjuiyje 
averaging over many microstructures are observed puu 
in particular, an increase in the average resistance at 
zero magnetic field called the weak-localization correc- 
tion (WLC).eI The two main features of these interfer- 
ence effects are their shape — the characteristic field or 
energy to which they are sensitive, or more precisely the 
functional dependence on these parameters — and their 
magnitude — simply how big these quantum corrections 
are compared to the classical conductance. Here we con- 
centrate on the magnitude and merely note thai ^thmiy 



for the shape has been extensively developedJ13'tlHl3li 

A theory of the magnitude of quantum interference ef- 
fects in chaotic cavities has recently been developed by 
making a statistical—atijiatz for the S'-matrix describing 
the scattering.BBllil'BEl Refs. ^ and |l6| developed a 



random S'-matrix theory by assigning to S an "equal 
a priori distribution" once the symmetry requirements 
were imposed. The results for the average, variance, 
and probability density of the conductance were in good 
agreement with numerical calculations for a chaotic cav- 
ity connected to two waveguides. M However, the random- 
matrix predictions for both the weak-localization correc- 
tion andp-the variance are larger than the experimental 
results.la'Ll In addition, the measured probability density 
is close to a Gaussian distributioncl when there are two 



propagating modes per lead (iV — 2), while random- 
matrix--theory predicts a Gaussian distribution only for 
iV>3.B 

Inherent in the model of Refs. Qq andHq is the assump- 
tion (among others) that one can neglect phase-breaking 
processes which destroy the coherence of the wave func- 
tion. In this paper we show that this assumption is 
largely responsible for the discrepancy between theory 
and experiment mentioned above. We make specific pre- 
dictions for the dependence of the quantum transport 
corrections on the degree of phase-breaking which may 
be tested by future experiments. 

To simulate the effects of phase-breaking events we 
adopt a model suggested by M. BiittikerO in addition to 
the physical leads 1, 2 attached to reservoirs at chemical 
potentials /ii, /i2, a- lead 3 connects the cavity to a phase- 
randomizing reservoir at /i3. This model has been dis- 
cussed extensively for disordered materialalj and, niore 
recently, for ballistic quantum dots by Marcus, et al.n A 
similar model has been used for abspsption of microwaves 
in chaotic-scattering from cavities.113 Requiring the cur- 
rent in lead 3 to vanish determines /ia; the two-terminal 
dimensionless conductance is then found to be 
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where Ty is the transmission coefficient for "spinless elec- 
trons" from lead j to lead i. The factor of 2 accounts for 
spin explicitly. We call TV the number of channels in leads 
1 and 2, N4, that in lead 3, and iVr = 2Ar -f N4,. 

We now make the fundamental assumption of an 
"equal a priori distribution" for the total Nt x Nt scat- 
tering matrix S, once the symmetry requirements have 
been imposed. S is, of course, unitary, and is symmet- 
ric in the absence of a magnetic field because of time- 
reversal symmetry. We assume that the statistics of the 
total S'-matrix are givpi by the Circular Ensembles of 
random matrix theoryE3 For B = the orthogonal en- 
semble (COE, denoted (3 — 1) is appropriate while for 



nonzero B we use the unitary ensemble (CUE, Af^ 2). 
In contrast to prewpus studies of the eigenphases,E3E3 tp-. 
tal transmissionjla'ta or individual S'-matrix elenients,t3 
we treat the statistics of g given in Eq. (|l]) . In the rest 
of the paper, we first derive results valid in the weak and 
strong phase-breaking limits, then combine these into an 
interpolation formula which simulations show to be valid, 
and finally compare with experiments. 

We start by recalling the result for the WLC and vari- 
ance at N^ = given in Ref. 03: 



dg = {g^^^^^ - (.g)('^='' = -iV/(27V + 1) 



AN{N +1) 



vBxg 



(2N+l)^{2N + 3) ' 

Ar2 



47V2_ I 



(3 = 1 
(3 = 2 



(2a) 



(2b) 



In addition to these results, the probability density of g 
was calculated for A^ = 1 — 3. 

The case N^ = 1 and /3 = 1 is a special one which 
can be analyzed in detail. From the joint jwobability 
distribution of S'-matrix elements in one rowJHil 

Nt 
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one can show that the joint distribution of the Tsj's is 

p.r^ .^ Ni2N-mn,n,)^-^sii-j:n,) ,.. 
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Remarkably, the {g) that one obtains by integrating over 
this distribution is identical to that for TV^ = 0. 

In the limit TV^ ^ 1, one can obtain a number of results 
to leading order in l/iV^. First, expand the conductance 
of Eq. (|i|) in powers of STij where Tij = (Tij) + STi,. By 
definition (STij) = 0; using the methods of Ref. E2, one 
finds that the correlations among the T^ 's are given by 

{ST,,6Tki) = Ap {N,N,{Nt + 5pi){NT + 25pi)5,k5ji 
- N,N,Ni{Nt + <5/5i)<5,fc - N.N.NkiNr + S„i)6,i 
+ N,N,NkNi{l + 6pi)} , (5a) 
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for i ^ j, k ^ I. (For /3 = 1, the indices of the T^ 's have 
been permuted so as to maximize coincidences.) Note 
that {STijSTki) is at least of order l/Nj. Thus 
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Here we haveta (Ty) = N,Nj/{Nt -f ^/3i), a result 
closely Eelated to the Hauser-Feshbach formula in nuclear 
physics.c3 The WLC is then given by 



5g = -N/N^ + 0{l/Nl). 



For the variance one finds from Eq, 
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and, for the ratio of the variances for (3 = 1,2, 
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Of course the magnitude of the quantum corrections de- 
creases as the number of phase-breaking channels in- 
creases. Note that the power controlling the decay of the 
variance is twice that of the average conductance. The 
deviation of the ratio of the variances from 2 is highly 
unusual; in fact, for A^ = 1 the ratio can be as high as 3. 
We present in Figs. 1, 2 and 3 results of simulations of 
the random-matrix model in order to check the asymp- 
totic behavior and suggest improvements. The points are 
obtained by generating random Nx x N^ unitary or or- 
thogonal matrices and computing g from Eq. ([^) . In Fig. 
1, the value N = 2 was selected for comparison with the 
experiment of Ref. 0. The log-log insets in Fig. 1 show 
that the convergence of the numerical to the asymptotic 
results [Eqs. (R)-(p|), dotted lines] is rather slow. The 
curves in the main parts of Fig. 1 are interpolation for- 
mulae suggested by the results for N^, = and iV^ :^ 1. 
For the WLC, Eqs. ( pa| ) and (Q) suggest the relation 
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For the fluctuations of the conductance, we combine 
square roots of variances. 
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because in the numerical simulations the change in vavg 
is clearly linear in N^, not quadratic. So, for (3 = 2 
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For (3 = 1 a, similar, but more complicated, expression 
holds. These interpolation formulae agree very well with 
the numerical results; ; the only significant deviation is 
for A^ = 1 and small N^ [Fig. 1(c)]. 

The results at fixed values of N^ shown in Fig. 2 are 
relevant to experiments at fixed temperature in which 
the size of the opening to the cavity is varied. Though 
6g and vavg are nearly independent of N in the per- 
fectly coherent limit — a result known as "universality" — 
phase-breaking channels cause the magnitudes to vary. 
Thus the universality can only be seen if N^ <ti N; oth- 
erwise, the behavior is approximately linear, as in some 



experimentsfl'tl Clearly, phase-breaking must be included 
in interpreting the experiments. 

In addition to the mean and variance, the probabil- 
ity density of the conductance, w{g), is experimentally 
measurable.u The numerical results in Fig. 3 show that 
as N^ increases w{g) tends towards a Gaussian and is 
therefore fully characterized by the mean and variance 
given abojse. For iV > 3 the distribution is essentially 
Gaussianta even for N^ = 0, so that N^ affects w{g), 
apart from changing the mean and variance, only for , 
3. For A^ = 1, the highly non-Gaussian distributions 
at N^ — become approximately Gaussian at N^ — 3; 
the intermediate distributions, N^f, = 1,2, are shown in 
Fig. 3. For iV = 2, the deviation of w{g) for N^ = 
from a Gaussian is smaller than for N = 1, and hence 
fewer phase-breaking channels produce a Gaussian dis- 
tribution. We show numerical results only for N^ = 1. 

Finally, we use our model to interpret the experiments 
of Refs. ^, ^ and y. First, the solid circles in Fig. 1 
show the measured N^, Sg, and varg of Ref. in the 
N = 2 case. Before comparison the experimental vari- 
ance must be corrected for thermal averaging: convolu- 
tion over the derivative of the Fermi function produces a 
reduction factor of ~ 0.22 — 0.38 for an electron temper- 
ature of 50 — 100 mK. We have increased the measured 
varg by the inverse of this reduction factor; the WLC is 
not affected by thermal averaging since it is already an 
average effect. The error bars shown result from both the 
uncertainty in temperature in the case of the variance and 
the experimental fluctuations at small B for the variance 
and WLC. (For the moment, we do not assign further 
physical significance to the oscillations seen in the exper- 
iment.) Note that in Fig. 1 we have not fit the theory to 
the data and yet the agreement is very good. Second, in 
Fig. 2 we show as solid symbols the data of all three Refs. 
^, g and ||. The value of A'^^ deduced from comparing to 
our calculations is in good agreement with the value esti- 
mated independently in the experiments: iV^ = 7 — 9 for 
Ref. I, iV0 = 2 for Ref. §, and iV^ = 4 - 8 for Ref. |. Fi- 
nally, the probability density of the conductance is found 
to be Gaussian in Ref. ^. For the estimated N^ = 4 — 8, 
we obtain a Gaussian distribution and so are consistent. 
Observation of the interesting non-Gaussian distribution 
of g obtained in the theory for N = 1,2 and N^ = 
requires greatly reduced phase-breaking. 

In conclusion, we have presented a random-matrix 
model that simulates the effects of phase-breaking events 
on electronic transport through ballistic chaotic cavities. 
The analysis of recent experiments indicates that one can 
find a value of iV^, the number of phase- breaking chan- 
nels, that allows a consistent description of the data. Fur- 
ther experiments are needed to test quantitatively the 
dependence on N and N^ that we predict. 

We thank CM. Marcus for several valuable conversa- 
tions. 

Note added — While preparing this paper for publica- 
tion, we received a preprint by P. W. Brouwer and C. W. 
J. Beenakker with some overlapping material for iV = 1. 
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FIG. 1. Magnitude of quantum transport effects as a 
function of the number of phase- breaking channels, A^^, on 
linear (main panels) and log-log (insets) scales, (a) The 
weak- localization correction (A'' = 2). (b) The variance in 
both the orthogonal (squares) and unitary (triangles) cases 
{N = 2). (c) The ratio of the variance in the orthogonal case 
to that in the unitary for A'^ = 1, 2, and 6; the arrows mark 
the N,/, — > oo limit for A'^ = 1, 2. Open symbols are numerical 
results (20, 000 matrices used, statistical error the size of the 
symbol). Solid lines are interpolation formulae. Dotted lines 
are asymptotic results. Solid circles are experimental results 
of Ref. H corrected for thermal averaging. The interpolation 
formulae are excellent except for A'^ = 1 and small A'^ [panel 



FIG. 2. The magnitude of quantum transport effects as a 
function of the number of channels in the leads. A', for fixed 
A'^^ = 0, 2, 4, and 8. (a) The weak-localization correction, 
(b) The variance for the orthogonal case {B = 0). (c) The 
variance for the unitary case (nonzero B). Open symbols are 
numerical results (as in Fig. 1). Solid lines are interpolation 
formulae. Solid symbols are experimental results of Refs. H 
(triangles) , M (squares) , and M (circles) corrected for thermal 
averaging. The introduction of phase-breaking decreases the 
"universality" of the results but leads to good agreement with 
experiment. 



FIG. 3. The probability density of the conductance in the 
orthogonal (first column) and unitary (second column) cases 
for A'^ = 1 (first row) and N = 2 (second row) . Increasing the 
phase-breaking from zero (dashed lines, analytic) to A^^ = 1 
(plus symbols, numerical) to A'^ = 2 (squares, numerical) 
moves the distribution towards a Gaussian. 
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Figure 3. 



